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This  paper  addresses  the  problem  of  controlling  the  receiver  aircraft  to  achieve  a  suc¬ 
cessful  aerial  refueling.  For  the  performance  verification  of  the  controller,  a  new  set  of 
nonlinear,  6-DOF,  rigid  body  equations  of  motion  for  the  receiver  aircraft  has  been  de¬ 
rived.  The  equations  are  developed  using  the  reference  frame  as  one  that  is  attached 
to,  translates  and  rotates  with  the  tanker  aircraft.  Furthermore,  the  nonlinear  equations 
contain  the  wind  effect  terms  and  their  time  derivatives  to  represent  the  aerodynamic  cou¬ 
pling  involved  between  the  two  aircraft.  These  wind  terms  are  obtained  using  an  averaging 
technique  that  computes  the  effective  induced  wind  components  and  wind  gradients  in  the 
receiver  aircraft’s  body  frame.  Dynamics  of  the  engine  and  the  actuators  are  also  included 
in  the  study.  A  linear  position-tracking  controller  has  been  designed  using  a  combination 
of  integral  control  and  optimal  LQR  design.  The  controller  does  not  use  the  information  of 
the  tanker’s  vortex  induced  wind  effects  acting  on  the  receiver  aircraft  as  well  as  the  mass 
change  that  occurs  during  refueling.  The  performance  of  the  controller  is  evaluated  in  the 
high  fidelity  simulation  environment  employing  the  new  sets  of  equations  of  motion.  The 
simulation  and  control  design  are  applied  to  a  tailless  fighter  aircraft  with  innovative  control 
effectors  and  thrust  vectoring  capability.  Various  allocation  schemes  for  redundant  control 
variables  are  analyzed  in  a  realistic  approach  maneuver  to  the  refueling  contact  position 
behind  the  tanker  aircraft.  In  this  paper,  the  performance  evaluation  is  presented  only 
during  the  initial  phase  of  aerial  refueling  maneuver  when  the  receiver  aircraft  maneuvers 
to  reach  the  refueling  contact  position. 

I.  Introduction 

Undoubtedly  the  most  demanding  and  critical  aircraft  missions  are  those  involving  multiple  vehicles  and 
especially  when  they  need  to  operate  in  close  proximity  with  each  other  such  as  in  formation  flight  and  aerial 
refueling  operations.  Two  of  the  most  significant  factors  that  affect  the  receiver’s  dynamics  in  an  aerial 
refueling  operation  are  the  time- varying  mass  and  inertia  properties  during  fuel  transfer,  and  the  wind  effect 
due  to  the  tanker’s  trailing  vortices.  While  the  effect  of  the  leader’s  vortices  on  the  follower  in  a  formation 
flight  is  beneficial  in  reducing  drag/fuel  consumed,* 1  the  influence  of  a  tanker’s  vortices  on  a  receiver  aircraft 
during  aerial  refueling2,3  can  be  detrimental  to  the  stability  of  the  receiver  aircraft. 

There  has  been  some  recent  work4-18  dealing  with  demonstration  of  the  benefit  of  and  issues  with  the 
control  system  development  for  aerial  refueling.  However,  they  do  not  study  the  effect  of  the  tanker’s  trailing 
wake  vortex  on  the  receiver  aircraft.  Most  of  the  previous  work  treat  the  vortex  as  an  unknown  disturbance 
or  stochastic  turbulence  in  the  control  law  design  procedure,  and  they  do  not  attempt  to  verify  their  control 
design  in  an  environment  where  the  aircraft  can  be  exposed  to  a  realistic  wake  vortex.  Since  vortex-induced 
velocities  acting  on  an  aircraft  are  highly  non-uniform,  standard  aerodynamic  force  and  moment  equations, 
based  on  airspeed,  angles  of  attack  and  sideslip,  and  uniform  wind  components  and  gradients  acting  at  the 
center  of  mass  (CM)  of  the  aircraft,  cannot  be  used  directly.  To  overcome  this  difficulty,  there  have  been 
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two  extreme  approaches  used  thus  far. 19-21  The  first  one  is  to  generate  a  database  of  induced  forces  and 
moments  for  a  specific  pair  of  aircraft  by  very  complicated  CFD  (Computational  Fluid  Dynamics)  models 
or  wind  tunnel  and/or  flight  test  measurements.  Obviously,  the  databases  are  extremely  costly  to  generate, 
too  computationally  intensive  to  use  for  even  near  real-time  simulation  and  very  difficult  to  use  in  control 
system  development.  Furthermore,  they  are  specific  to  certain  pairs  of  aircraft  and  flight  condition.  Even 
in  this  case,  incorporation  of  the  nonuniform  wind  effect  in  the  equations  of  motion  as  induced  force  and 
moment  terms  invalidates  the  definition  of  aero-related  variables.  For  example,  airspeed  is  the  speed  of  the 
vehicle  relative  to  the  air  and  it  is  directly  affected  by  the  presence  of  wind,  which  is  ignored  when  the  effect 
of  wind  is  only  included  as  an  additional  force  or  moment.  The  other  extreme  is  either  to  completely  avoid 
modeling  the  effect  of  nonuniform  wind  and  model  it  as  an  unknown  disturbance  or  to  model  it  in  such  a 
simplified  manner  that  it  is  not  reliable  to  use  in  verification  work. 

Recently,  the  authors  have  developed  a  novel  modeling  technique  that  enables  the  use  of  standard  dynamic 
equations  of  motion  and  aerodynamic  build-up  equations  with  wind  effect  terms  included.3,22-28  First 
step  is  to  derive  the  equations  of  motion  of  an  aircraft  in  which  the  wind  effect  is  explicitly  stated  under 
the  assumption  that  the  aircraft  is  is  a  uniform  wind  field,  i.e.  uniform  wind  components  and  uniform 
wind  gradients.  Since  the  vortex-induced  wind  field  acting  on  the  encountering  aircraft  is  non-uniform 
in  nature,  the  novel  idea  is  to  approximate  the  non-uniform  induced  wind  components  and  gradients  by 
equivalent  uniform  wind  and  gradients.  Once  a  reasonable  approximation  is  achieved,  the  implementation  of 
aerodynamic  coupling  between  aircraft  flying  in  close  proximity  becomes  far  more  direct  and  computationally 
efficient  than  the  conventional  procedure  since  the  necessity  for  calculating  the  induced  forces  and  moments 
during  the  simulation  is  overcome.  This  approach  has  been  proven  very  useful  and  accurate  in  the  case  of 
formation  flight  modeling  and  control.22-27 

In  this  paper,  the  same  vortex  effect  modeling  technique  is  applied  to  the  aerial  refueling  problem.  First, 
the  equations  of  motion  of  the  receiver  aircraft  are  derived  in  terms  of  the  position  and  orientation  of  the 
receiver  relative  to  the  tanker  aircraft  while  the  velocity  vector  of  the  receiver  is  represented  as  the  sum  of 
the  velocity  of  the  aircraft  relative  to  the  surrounding  air  and  the  velocity  of  the  air  relative  to  an  inertial 
frame.  Then,  the  full  nonlinear  equations  of  motion  including  the  terms  with  translational  and  angular 
wind  velocity  components  are  implemented  in  a  MATLAB/Simulink  environment  along  with  a  vortex  model 
employing  the  vortex  effect  modeling  technique.  The  simulation  model  is  developed  specifically  for  a  pair 
of  tanker  receiver  aircraft:  KC-135R  as  the  tanker  and  a  tailless  fighter  aircraft  with  innovative  control 
effectors  (ICE)  and  thrust  vectoring. 

For  the  receiver  to  track  a  commanded  trajectory  to  approach  the  refueling  contact  position  behind  the 
tanker,  a  multi-input  multi-output  (MIMO)  linear  controller  is  designed  and  implemented  in  the  simulation 
environment.  The  performance  of  the  controller  is  evaluated  in  the  simulation  environment  with  the  full 
nonlinear  equations  and  trailing  wake  vortex  from  the  tanker  aircraft.  Evaluation  of  the  controller  is  carried 
out  in  three  different  allocation  schemes  of  the  redundant  control  variables:  three  control  effectors,  thrust 
magnitude  and  two  thrust  vectoring  angles. 

II.  Modeling  of  the  Receiver  Dynamics  Relative  to  the  Tanker 

In  an  efficient  aerial  refueling  operation,  the  receiver  aircraft  needs  to  be  controlled  with  respect  to  the 
tanker’s  position  and  orientation  rather  than  with  respect  to  an  inertial  reference.  Additionally,  the  receiver 
aircraft  will  be  exposed  to  a  nonuniform  wind  field  during  the  whole  refueling  operation  when  it  is  in  the 
proximity  of  the  tanker  due  to  the  trailing  vortex  of  the  tanker.  Thus,  a  new  set  of  nonlinear  equations  are 
derived3  and  used  in  this  study  to  represent  the  position  and  orientation  of  the  receiver  relative  to  the  tanker 
and  at  the  same  time  to  explicitly  represent  the  vortex  effect  on  the  dynamics  of  the  receiver  aircraft. 

A.  Translational  Kinematics  and  Dynamics 

Since  we  are  interested  in  the  position  of  the  receiver  relative  to  the  tanker  and  not  its  absolute  position 
with  respect  to  the  ground,  the  translational  kinematics  equation  is  written  in  terms  of  the  position  vector 
of  the  receiver  with  respect  to  tanker  (Fig.  1 )  : 

C  R-BftWFt  f/j  T  If  RbtI  ^ Bt  T  ^(^BjrBt  T  R-BuBt^Bt  )C  (1) 

where  £  is  the  position  of  the  receiver  relative  to  the  tanker  expressed  in  the  body  frame  of  the  tanker, 
Rbrwr  is  the  rotation  matrix  from  the  receiver  wind  frame  to  body  frame,  Vw  is  the  velocity  of  the  receiver 
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Figure  1.  Intermediate  reference  frames:  Tanker  and  Receiver’s  body  frame 


relative  to  the  surrounding  air  expressed  in  the  receiver  wind  frame,  W  is  the  velocity  of  the  surrounding  air 
relative  to  the  ground  expressed  in  the  receiver  body  frame,  Rbrbt  is  the  rotation  matrix  from  tanker  body 
frame  to  receiver  body  frame,  Rbti  is  the  rotation  matrix  from  the  inertial  frame  to  tanker  body  frame  and 
Tbt  is  the  velocity  of  the  tanker  relative  to  the  inertial  frame.  Also,  note  that  S(-)  is  the  skew-symmetric 
operation  on  the  representation  of  a  vector  and  defined  as 


s0) 


0 

x3 

-X2 

-X3 

0 

Xi 

X2 

-xi 

0 

(2) 


for  an  arbitrary  vector  x  with  the  representation  [xi  X2  x^]7* . 

The  equation  of  translational  dynamics  of  the  receiver  aircraft  including  the  wind  effect,  in  matrix  form, 
is  given  below: 


’  Vr  ' 

pR 

—  °n 

<3r 

^(^BrBt)  T  ^BrBt  ^(^Bt  )^BrBt 


R-BrWr.  Vw  +  W 


H - -fA1 1  RbrbtR-bti-Mi?  +  Rbrw»  Ar  +  PR 
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where 


C-  1 

--R 


cos  a  cos  /?  sin  (3 

—  cos  a  sin  (3  pk  cos  f3 

—  p-  sec  (3  sin  a  0 


cos  / 3  sin  a 
—  pk  sin  a  sin  (3 
p-  cos  a  sec  j3 


(4) 


In  this  form  of  the  translational  equations,  all  the  translational  and  rotational  states  of  the  receiver  are 
with  respect  to  the  tanker  body  frame.  Furthermore,  the  position  vector  is  represented  in  the  tanker  body 
frame.  Another  advantage  of  this  representation  is  due  to  the  fact  that  the  motion  of  the  tanker  aircraft 
-  both  translational  and  rotational  -  are  passed  as  exogenous  inputs  to  the  receiver  aircraft.  The  variables 
included  in  this  category  are  velocity  ,  rgT ,  orientation  in  terms  of  euler  angles  i/t,  </>t ,  and  angular 

velocity  pt ,  9t,  rr  and  translational  velocity  of  the  tanker,  all  relative  to  the  inertial  frame. 
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The  external  forces  acting  on  the  receiver  are  the  gravitational  force  vector  (expressed  in  the  inertial 
frame),  the  aerodynamic  force  vector  AR  (expressed  in  the  wind  frame  of  the  receiver)  and  the  propulsive 
force  vector  PR  (expressed  in  the  body  frame  of  the  receiver).  In  general,  the  representations  of  the  force 
vectors  M R,  AR  and  PR  are 


0 

-Dr  ' 

’  Tx  ' 

Mr  = 

0 

Ar  = 

-Sr 

Pr  = 

Ty 

.  m9  _ 

—  Lr 

Tz 

The  propulsion  force  -PR  has  three  components,  which  are  functions  of  thrust  magnitude  TR  and  the  orien¬ 
tation  of  the  thrust  vector.  The  thrust  vectoring  is  parameterized  by  the  angle  of  the  thrust  vector  with  the 
receiver’s  body  xy-  and  xz-  planes.  Thus,  as  seen  in  Fig.  21,  the  components  of  the  thrust  are 

Tx  =  Tr  cos  Sz  cos  Sy 

Ty  =  Tr  sin  Sz  (6) 

Tz  =  Tr  cos  5z  sin  Sy 

Note  a  positive  Sy  rotation  of  the  thrust  generates  a  positive  thrust  component  in  the  positive  ^-direction 
while  inducing  a  positive  pitching  moment.  Similarly,  a  positive  Sz  rotation  of  the  thrust  generates  a  positive 
thrust  component  in  the  positive  y-direction  while  inducing  a  negative  yawing  moment. 


Figure  2.  Thrust  vectoring  angles  and  moment  arms 
The  aerodynamic  forces  are  given  by 

Dr  =  \PV2SRCD,  (7) 

Sr  =  \pV2SRCs,  (8) 

Lr  =  \pV£SRCL,  (9) 

where  SR  is  the  reference  area  of  the  receiver  and  p  is  the  ambient  air  density.  The  aerodynamic  coefficients 

are 

Cd  =  Cdo  +  Cdoc®  +  CDa20t2  +  CDSe5e  +  CD6e25l  (10) 

Cs  =  Cso  +  CspP  +  Cs8aSa  +  CssrSr  (11) 

Cl  =  Clo  +  CLa&  + CLa.2i.0t  —  aref)2  +  C'Lq-^Tj-qrel  +  CL8c5e  (12) 


where  a  is  the  angle-of-attack,  /?  is  the  side  slip  angle,  qrei  is  the  angular  velocity  of  the  receiver  relative 
to  the  surrounding  air  around  the  body-y  axis  and  (Sa,  Se,  Sr )  are  the  deflections  of  the  control  effectors 
(aileron,  elevator,  rudder)  as  the  conventional  control  surfaces  or  (elevon,  pitch  flap,  clamshell)  as  in  the  ICE 
(Innovative  Control  Effectors)  aircraft,  respectively. 
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B.  Rotational  Kinematics  and  Dynamics 

The  rotational  motion  of  the  receiver  aircraft,  similar  to  its  translational  motion,  is  also  analyzed  with 
reference  to  the  tanker  body  frame.  Even  though  the  standard  rotational  kinematics  and  dynamics  equations 
are  used,  their  interpretations  are  different  because  both  the  angular  position  and  velocity  of  the  receiver 
aircraft  are  relative  to  the  tanker  body  frame,  an  accelerating  and  rotating  reference  frame. 

The  rotational  kinematics  equation  in  the  matrix  form  is  the  well  known  standard  equation: 

S(ubrbt)  ;  (1^) 

where  u>brbt  is  the  representation  of  the  angular  velocity  vector  of  the  receiver  aircraft  relative  to  the  tanker 
body  frame  expressed  in  its  own  body  frame  as 


^  B  rBt 


PRt 

Q.Rt 

rRT 


(14) 


The  matrix  form  of  the  rotational  dynamics  of  the  receiver  is  also  modeled  with  the  standard  rotational 
dynamics  equation: 


^BRBrp  Xpp  ~ k  Xpp  ®(^Bj^Bx  ” k  ^^B ^Bj ^B^ )— R,  — k  (1A) 

where  IR  is  the  inertia  matrix  of  the  receiver  aircraft,  Mbr  is  the  moment  of  the  external  forces  around  the 
origin  of  the  receiver  body  frame  and  expressed  in  the  receiver  body  frame  as 


Mb 


R 


c 

M 


(16) 


.  The  moment  has  two  main  components;  due  to  aerodynamic  forces  and  due  to  the  thrust,  thus 


£  = 

-pVftSRbC£  —  £±ZTy  +  AyTz 

(17) 

M  = 

-pVrSrcCm  ~  AzTx  -  A XTZ 

(18) 

M  = 

2  pCrS  r})C m  —  AyTx  +  A  xTy 

(19) 

where  b  is  the  wingspan,  c  is  the  cord  length  of  the  receiver  aircraft,  and  (A^,,  Aj,,  Az)  are  the  moment  arms 
of  the  thrust  in  the  body  frame  of  the  receiver  (see  Fig.  2).  The  aerodynamic  moment  coefficients  are 


Cc 

CM 

Cm 


Cco  +  C £8  a  5  a  +  C £8  r  Sr  +  C£pf3  +  C£p  Prel  +  C £r——rrel 

IV R  IV  R 

Q 

CmO  +  C£aa  +  C£SeSe  +  C'MqR^qrel 

Z VR 

Cm  0  +  Cm  5  a  &a  +  Cm&t^t  +  CmpP  +  CVp  prel  +  CMrMrrTrel 

ZVr  ZVr 


(20) 

(21) 

(22) 


where  (preBQreiii~rei)  are  components  of  the  angular  velocity  of  the  aircraft  relative  to  the  surrounding  air. 
When  the  aircraft  is  in  a  vortex  field  as  in  the  case  of  tanker’s  trailing  wake  vortex  field,  these  angular  velocity 
components  will  be  different  from  the  angular  velocity  relative  to  the  tanker  of  relative  to  the  inertial  frame. 
Their  calculation  will  be  explained  in  Section  III.  There  will  be  some  additional  dynamics  induced  on  the 
receiver  due  to  the  facts  that  the  actual  fuel  transfer  will  in  fact  transfer  momentum  and  that  the  inertia 
properties  of  the  receiver  aircraft  will  be  time-varying.  These  additional  dynamic  effects  both  translational 
and  rotational  -  are  not  considered  in  the  present  study. 
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C.  Engine  Dynamics 

The  thrust  generated  by  the  engine  (T)  is 

T  =  £  Tmax  (23) 

where  £  denotes  the  instantaneous  throttle  setting  and  Tmax  is  the  maximum  available  thrust  and  assumed 
to  be  constant  in  this  paper.  The  engine  dynamics  is  modeled  as  that  of  a  first  order  system  with  time 
constant  r.  Therefore,  we  have  .  F  —  Ft 

Z  =  ,  (24) 

T 

where  Zt  is  the  commanded  throttle  setting  (0<  Zt  <  1). 

D.  Actuator  Dynamics 

For  the  present  study,  only  the  actuator  saturation  and  rate  limit  effects  are  considered.  Other  dynamics 
would  be  included  in  the  future  work.  The  deflection  range  attainable  from  the  elevon  is  (-30  deg,  30  deg), 
from  the  pitch  flap  (-30  deg,  30  deg)  and  from  the  clamshells  (0  deg,  60  deg).  All  the  three  control  effectors 
have  a  rate  limit  of  ±90  deg/sec.  Likewise,  the  thrust  vectoring  has  a  limit  of  ±  30  deg  in  both  directions 
and  a  rate  limit  of  ±30  deg/sec. 

III.  Modeling  the  vortex  and  its  effect 

It  is  to  be  noted  that  the  wind  effect  terms  constituting  the  elements  W,  W  in  the  receiver’s  equations  of 
motion  presented  earlier  are  considered  to  be  based  on  the  uniform  wind  distribution  acting  at  the  receiver’s 
CM,  expressed  in  its  body  frame.  But,  the  vortex-induced  wind  field  acting  on  the  receiver  aircraft  is  non- 
uniform  in  nature.  Therefore,  to  be  able  to  use  the  standard  aircraft  equations  of  motion  without  doing 
any  modifications,  there  is  a  need  to  approximate  the  non-uniform  induced  wind  components  and  gradients 
by  equivalent  uniform  wind  and  gradients.  Once  a  fairly  reasonable  approximation  can  be  achieved,  the 
implementation  of  aerodynamic  coupling  between  the  tanker  and  the  receiver  becomes  far  more  direct  and 
computationally  efficient  than  the  conventional  procedure  which  involves  first  the  calculation  of  induced 
forces  and  moments  from  the  wind  distribution,  and  then  inserting  these  forces  and  moments  in  the  aircraft 
dynamics  equations.  The  need  for  a  simple  and  fairly  accurate  method  of  approximating  the  non-uniform 
vortex-induced  wind  field  by  its  uniform  equivalent  forms  the  motivation  for  the  material  presented  in  this 
section. 


Figure  3.  Trailing  vortex  from  the  wings  and  horizontal  tail 

In  our  dynamic  model  for  aerial  refueling,  the  tanker  is  considered  to  produce  two  pairs  of  straight,  semi- 
infinite  trailing  vortex  filaments  -  one  from  the  wings  and  one  from  the  horizontal  tail  -  that  induce  additional 
wind  velocities  on  the  body  of  the  receiver  aircraft  (see  Fig.  3).  These  vortex-induced  wind  velocities  cause 
changes  in  the  forces  and  moments  experienced  by  the  receiver.  However,  instead  of  attempting  to  directly 
estimate  the  induced  forces  and  moments  on  the  follower,  the  induced  wind  velocities  and  wind  gradients 
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are  computed.  The  induced  wind  velocities  are  written  as  a  function  of  the  relative  separation  as  well  as  the 
relative  orientation  between  the  tanker  and  the  receiver  using  a  modified  horseshoe  vortex  model  based  on  the 
Helmholtz  profile.  Since  the  induced  wind  and  wind  gradients  are  non-uniform  along  the  body  dimensions  of 
the  receiver  aircraft,  an  averaging  technique  is  implemented  to  compute  the  effective  wind  and  wind  gradient 
as  uniform  approximations.  The  effective  wind  components  and  gradients  are  introduced  into  the  nonlinear 
aircraft  equations  that  include  the  components  of  wind  and  the  temporal  variation  of  wind  in  the  body  frame 
to  determine  the  effect  on  the  receiver’s  dynamics.  The  effect  of  vortex  decay  over  time  is  also  included  in 
our  model.  Special  care  has  been  taken  to  accommodate  different  geometrical  dimensions  for  the  tanker  and 
the  receiver  aircraft  and  also  to  include  many  useful  geometrical  parameters  of  the  aircraft  like  the  wing 
sweep  angle,  the  dihedral  angle  and  the  relative  distance  between  the  center  of  mass  of  the  aircraft  and 
the  aerodynamic  center  of  the  wing,  in  estimating  the  vortex-effect  experienced  during  aerial  refueling.  The 
previous  publication3, 22-28  of  the  authors  should  be  referred  to  for  further  details  of  the  actual  vortex  model 
and  the  averaging  technique  used  to  estimate  the  vortex-effect  on  the  receiver. 

IV.  Control  Design 

Once  the  approach  and  fly-away  trajectories  are  determined,  a  position-tracking  controller  is  needed  to 
make  the  receiver  follow  the  reference  trajectories.  During  the  station-keeping  phase  of  the  refueling  flight, 
the  reference  trajectory  of  the  receiver  is  a  constant  with  respect  to  the  tanker’s  body  frame.  In  the  approach 
and  fly-away  maneuvers,  the  reference  trajectories  should  be  generated  to  ensure  the  overall  safety  of  the 
receiver  and  the  tanker.  Therefore,  it  is  very  important  to  have  a  controller  that  can  fly  the  receiver  close 
to  the  reference  trajectories. 

The  primary  requirement  of  the  control  design  is  the  tracking  of  the  generated  trajectories,  with  zero 
steady-state  error  in  the  x,  y,  z  coordinates  in  the  tanker’s  body  frame,  under  the  disturbance  of  trailing 
vortex,  time  variation  of  the  inertia  properties  of  the  receiver  and  the  possible  steady  maneuvers  of  the 
tanker’s  body  frame.  Meanwhile,  the  control  inputs  generated  by  the  controller  should  not  cause  significant 
saturation  on  the  magnitudes  and  rates  of  the  actuators.  Moreover,  during  the  transient,  overshoot  or 
undershoot  on  trajectory  response  should  be  minimized  to  ensure  the  safety  of  the  refueling.  At  the  same 
time,  the  response  of  the  closed  loop  system  should  be  fast  enough  so  that  the  approach  and  fly-away 
maneuvers  are  completed  as  planned  and  the  high-wind  regions  of  the  trailing  vortex  field  are  exited  in  a 
timely  fashion.  Additionally,  during  the  approach,  fly-away  and  station-keeping  maneuvers,  the  angle-of- 
attack  and  the  airspeed  should  not  be  close  to  their  corresponding  stall  values.  In  this  regard,  very  big  pitch 
angle  should  not  be  commanded.  Finally,  to  ensure  the  safety  of  the  aircraft,  the  bank  angle  should  be  small 
relative  to  its  nominal  value. 

A  combination  of  a  multi-input  multi-output  LQR  and  integral  control  is  employed  in  designing  the 
position  tracking  controller  satisfying  the  above  requirements.  First,  the  nominal  conditions  of  the  receiver 
aircraft  are  calculated  at  the  steady  level  flight  condition  that  the  tanker  is  flying,  i.e.  the  airspeed  is  equal 
to  the  tanker’s  speed,  zero  bank  and  zero  yaw  angles  (recall  that  the  euler  angles  are  relative  to  the  tanker) 
and  zero  side  slip  angle.  Then,  the  nonlinear  equations  of  motion  of  the  receiver  presented  in  Section  II 
are  linearized  around  this  nominal  condition.  Note  that  in  the  linearization  procedure,  the  wind  terms 
are  ignored  while  in  the  simulations  the  full  nonlinear  equations  of  motions  are  employed.  The  linearized 
equations  of  motion  of  the  receiver  in  state-space  form  is  obtained  as 

x  =  Ax  +  Bit  +  Hrc  (25) 

where  A  €  3?12x12,  B  G  3?12x6,  H  G  5ft12x6  (for  the  numerical  values  of  the  matrices,  see  Appendix  A),  the 
state  vector  is 

x  =  [AV  A/3  Aq  A p  A q  A r  Aip  A 9  Acj>  Ax  Ay  A z]T  (26) 

the  control  input  vector  is 

u  =  [Ai5a  A 5e  A Sr  A£  ASy  A<52]T  (27) 

and  the  disturbance  vector  due  to  the  motion  of  the  tanker  is 

u=  [AVxt  A Vyt  AVZ T  Ax[)t  A Ot  A (28) 

Since  the  position  tracking  controller  is  to  be  designed  for  the  receiver  relative  to  the  tanker,  the  outputs 
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to  be  tracked  are  (Ax,  Ay, Ax).  Thus,  output  vector  is  chosen  to  be 

y  =  [Ax  Ay  A z]T  (29) 

To  ensure  zero  tracking  error  at  steady  state  condition,  the  state  space  equations  are  augmented  by  three 
integrators,  one  for  each  position  error: 

e  =  y  -  Vc_  (30) 

where  yc  =  [Axc  A yc  Axc]T  is  the  commanded  trajectory  for  the  receiver  in  the  body  frame  of  the  tanker 
for  approaching  the  refueling  contact  position.  Thus,  in  scalar  form 


ex  =  Ax  —  Axc 

ey  =  Ax  —  Axc  (31) 

ex  =  Ax  —  Axc 


By  including  the  augmentation  states  in  the  state-space  equations,  the  augmented  state  equation  becomes 

(32) 
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Figure  4.  State  feedback  and  integral  control  structure 


Using 

function: 


LQR  design  technique,  the  state  feedback  gain  matrix  [Kx  Ke]  is  obtained  to  minimize  the  cost 


OO 


0 


where  Q  £  3ff15x15  is  symmetric  positive  semidefinite  and  R  £  3?6x6  is  symmetric  positive  definite.  Thus, 
the  state  feedback  with  the  integral  control  is 


u  =  —  Kx  x  —  Ke  e 


(34) 


V.  Simulation  Results 

The  Nonlinear  equations  of  motion  of  the  receiver  are  implemented  in  a  MATLAB/Simulink-based  model 
along  with  a  tanker  model.  Both  the  receiver  and  the  tanker  models  are  equipped  with  their  own  controllers. 
The  receiver  controller  is  designed  and  implemented  in  such  a  way  that  various  levels  of  control  input 
allocations  can  be  tested.  The  control  variable  allocation  is  achieved  by  varying  the  respective  elements  of 
the  control  weighting  matrix  R  of  the  cost  function  in  Eq.  33.  Three  different  cases  are  simulated  and 
compared  in  this  section: 

Case— 1  :  A  combination  of  control  effectors  and  thrust  vectoring  is  used.  R  is  chosen  to  a  diagonal  matrix 
with  the  following  elements 

R(m)  =  (1,1,10,1000,50, 10)  (35) 
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Case— 2  :  Only  control  effectors  are  used  without  any  thrust  vectoring.  R  is  chosen  to  a  diagonal  matrix 
with  the  following  elements 


R(i,  i)  =  (1, 1, 10, 1000, 1000000, 1000000) 


(36) 


Case— 3  :  Pitch  Flap  and  Clamshell  are  fixed  at  their  nominal  values.  Elevon  and  thrust  vectoring  are  only 
control  variables  used.  R  is  chosen  to  a  diagonal  matrix  with  the  following  elements 


R(i,  i)  =  (10, 1000000, 1000000, 1000, 50, 100) 


(37) 


In  all  the  three  cases,  Q  is  chosen  to  a  diagonal  matrix  with  the  following  elements 


Q(M)  =0.1  x  (1,1,1,1,1,1,1,1,1,1,0.001,1,1,0.1,1) 


(38) 


The  refueling  is  performed  at  altitude  of  7010  meters  while  the  tanker  is  flying  at  speed  of  200  ml  sec. 
The  simulations  start  when  the  receiver  is  laterally  offset  from  the  tanker  by  60.96  meters,  longitudinally 
15.24  meters  behind  the  refueling  contact  position  and  vertically  at  the  same  level  as  the  the  contact  position. 
The  refueling  position  is  25.33  meters  directly  behind  and  6.46  meters  below  the  c.g.  of  the  tanker  without 
any  lateral  offset.  In  other  words,  the  receiver  is  initially  at  (-40.56,60.96,6.46)  in  the  body-frame  of  the 
tanker  and  should  go  to  and  stay  at  (-25.33,0,6.46)  for  refueling  .  The  receiver  should  move  from  the  initial 
position  to  the  refueling  contact  position  by  first  maneuvering  laterally  right  behind  the  tanker  and  then 
moving  forward  to  the  contact  position  without  any  altitude  change  (see  Figs.  5  and  6). 

Figure  5  shows  the  x-  and  y-  components  of  both  the  commanded  and  actual  trajectory  of  the  receiver 
in  the  body  frame  of  the  tanker  in  case-1.  Fig.  6  illustrates  the  x-  ,  y-  and  z-  components  of  the  trajectory 
in  time  domain.  As  will  be  shown  more  clearly,  when  the  simulation  starts,  the  vortex-induced  wind  is  not 
on  and  at  10  seconds,  the  wind  is  turned  on  and  gradually  increased  to  the  normal  level.  This  is  done  to 
ensure  the  start  of  the  simulation  without  any  numerical  problem  and  also  to  see  the  effect  of  the  vortex  in 
the  initial  flight  configuration.  As  seen  Fig.  6,  the  lateral  maneuver  is  initiated  at  around  50  seconds  and 
lasts  about  a  minute  until  the  receiver  is  directly  behind  the  tanker  at  110  seconds.  Then,  at  125  seconds,  the 
receiver  starts  moving  forward  towards  and  reaches  the  the  refueling  contact  position  in  about  50  seconds. 
The  trajectories  of  the  other  two  cases  are  not  presented  because  they  all  are  similar  to  each  other.  Both  the 
figures  show  that  the  commanded  trajectories  are  closely  tracked  during  the  maneuver.  Note  that  positive 
z-direction  in  down.  Thus,  the  two  spikes  in  the  z-plot  in  Fig.  6  are  in  fact  due  to  slight  altitude  lost  first 
when  the  lateral  maneuver  starts  and  second  when  the  forward  motion  is  initiated. 


Commanded  and  Actual  Trajectory  in  xy-plane  in  case-1 


Figure  5.  Trajectory  of  the  receiver  relative  to  the  tanker  in  xy-plane 

Figures  7  8  and  9  shows  the  initial  levels  of  the  control  variables,  how  they  vary  during  the  entire  maneuver 
and  their  final  levels  at  the  refueling  positions  in  all  the  three  cases.  Fig.  7  illustrates  the  deflections  of  the 
three  control  effectors,  elevon,  pitch  flap  and  clamshell.  Note  that  the  pitch  flap  and  clamshell  stay  constant 
in  case-3  because  they  are  fixed  at  their  nominal  values.  In  all  the  three  cases,  very  small  deflections  and 
deflection  rates  (much  smaller  than  their  respective  saturation  levels)  in  all  the  three  effectors  are  used  to 
move  the  receiver  to  the  refueling  position.  Both  elevon  and  clamshell,  as  expected,  are  mainly  used  during 
the  lateral  maneuver  and  stays  almost  constant  during  the  forward  motion.  The  behavior  of  the  pitch  flap 
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Commanded  and  Actual  x,y  and  z  positions  [m]  in  case-1 


Figure  6.  Commanded  and  actual  x,  y  and  z  —  positions  of  the  receiver 


is  significantly  different  between  case-1  and  case-2.  In  case-2,  it  is  used  only  during  the  lateral  maneuver 
but  in  case-1  used  in  both  lateral  and  forward  maneuvers.  This  is  because  in  case-1,  it  is  coupled  with  the 
thrust  vectoring  as  can  be  seen  in  Fig.  8. 


C/3 

E 

J2 

O 


Deflection  of  Control  Effectors  [deg] 


Figure  7.  Deflections  of  the  three  control  effectors  in  the  three  cases 


Figure  8  shows  the  rotation  of  the  thrust  vector  in  terms  of  its  angles  with  the  xy-  and  xz-  planes  of 
the  body-frame  of  the  receiver.  Note  that  rotation  of  angle  with  xy-plane  induces  pitching  moment  while 
generating  vertical  force.  Similarly,  rotation  of  angle  with  xz-plane  induces  yawing  moment  while  generating 
lateral  force.  Also  note  that  both  the  angles  stay  zero  in  case-2  since  it  represents  no-thrust- vectoring 
configuration.  In  case-1,  there  is  a  direct  correlation  between  angle  with  xy-plane  and  pitch  flap  since  they 

10  of  15 


American  Institute  of  Aeronautics  and  Astronautics 


0  50  100  150  200 

time  [sec] 


Figure  8.  Thrust  vectoring  angle  variation  in  the  three  cases 


both  mainly  generate  pitching  moment.  In  case-3,  since  the  pitch  flap  is  fixed  and  thus  no  coupling,  the 
angle  with  xy-plane  increases  during  the  lateral  maneuver  and  stays  at  this  new  level  afterwards.  The  angle 
with  xz-plane  in  case-1  is  decreased  slightly  while  the  clamshell  angle  is  increased.  However,  in  case-3, 
it  is  increased  to  about  6  degrees  since  the  clamshell  is  not  used  at  all.  When  the  clamshell  deflections  in 
cases  1-  and  2-  are  compared,  it  can  be  seen  that  its  deflection  is  higher  in  case-1  and  because  it  needs  to 
compensates  the  moment  induced  by  the  rotation  of  the  thrust  by  the  negative  angle  with  xz-plane. 

Figure  9  illustrates  the  level  of  throttle  in  all  the  three  cases.  First  note  that  when  the  vortex  is  turned 
on  at  10  seconds ,  the  level  of  throttle  decreases  slightly  since  the  receiver  aircraft  is  in  the  upwash  region  of 
the  tanker’s  trailing  vortex.  During  the  maneuver  in  all  the  cases,  the  variation  of  the  throttle  setting  is  very 
similar.  Once  the  receiver  is  at  the  refueling  position,  the  throttle  level  in  all  the  cases  are  much  higher  than 
that  at  the  initial  position.  This  is  because  at  the  final  position,  the  aircraft  is  subject  to  strong  downwash, 
which  results  in  higher  required  thrust.  The  slight  difference  in  the  final  level  of  throttle  between  the  cases  is 
due  to  the  different  thrust  vectoring  angles  and  different  orientation  of  the  aircraft  at  the  refueling  position. 


Figure  9.  Variation  in  throttle  level  during  the  maneuver  in  the  three  cases 


Figure  10  shows  how  the  orientation  of  the  receiver  changes  during  the  maneuvers  in  all  the  three  cases 
in  terms  of  euler  angles  relative  to  the  tanker’s  body  frame.  Note  that  during  the  whole  maneuver,  all  the 
euler  angles  are  kept  small  by  the  controller  while  the  position  of  the  receiver  aircraft  changes  and  the  level 
of  vortex  exposure  increases.  However,  there  is  a  nonzero  yaw  angle  and  it  is  not  corrected  by  the  controller 
in  all  the  three  cases.  This  is  expected  because  the  controller  is  designed  to  track  the  position  with  zero 
steady-state  error  by  the  additional  integral  controllers  for  the  three  commanded  position  signals  while  the 
LQR  controller  only  guarantees  the  stability.  However,  three  more  integral  control  terms  could  be  added  to 

11  of  15 


American  Institute  of  Aeronautics  and  Astronautics 


Euler  Angles  [deg] 


/  |- . r  ■  ■  . , . | . — 

0  50  100  150  200 

time  [sec] 


Figure  10.  Orientation  of  the  aircraft  relative  to  the  tanker  in  the  three  cases 


the  controller  since  there  are  six  control  variables  in  this  receiver.  Note  that  the  final  lateral  equilibrium 
condition  involves  the  highest  yaw  angle  and  clamshell  deflection  in  case-1,  slightly  smaller  ones  in  case-2 
and  smallest  yaw  angle  while  biggest  thrust  angle  with  xz-plane  in  case-3.  On  the  other  hand,  there  is  a 
nonzero  steady-state  bank  angle  in  case-3. 


Induced  Wind  Components  [m/sec]  and  Gradients  [deg/sec] 
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Figure  11.  The  wind  components  and  gradients  the  receiver  is  exposed  to  in  case— 1 

The  vortex  induced  wind  components  and  gradients  are  shown  in  Fig.  11  in  case-1.  In  the  other  two 
cases  the  plots  are  similar.  First  note  that  at  10  seconds ,  the  vortex  is  ’’turned  on”  and  the  aircraft  is  exposed 
to  small  upwash  and  sidewash  and  a  slight  ’’rolling”  gradient.  During  the  lateral  maneuvers,  as  the  receiver 
gets  laterally  closer  to  the  tanker,  both  the  amounts  of  effective  wind  components  and  gradients  increases. 
At  about  70  seconds  when  the  lateral  distance  to  the  tanker  is  about  %60  percent  of  the  tanker’s  wingspan, 
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the  highest  upwash  is  experienced.  As  the  receiver  gets  closer  to  the  tanker  when  the  lateral  distance  to 
the  tanker  is  about  %45  of  the  tanker’s  wingspan,  the  receiver  experiences  the  greatest  rolling  vortex  and 
the  upwash  turns  into  downwash  while  the  sidewash  increases  dramatically.  This  is  manifested  in  the  rolling 
oscillation  as  seen  in  Fig.  10  and  altitude  drop  in  Fig.  6.  During  this  transition,  a  yawing  vortex  gradient 
is  experienced  while  a  pitching  gradient  develops.  As  the  receiver  gets  even  closer  to  a  position  right  behind 
the  tanker,  the  rolling  vortex  gradient  and  sidewash  disappear  and  downwash  increases  to  its  highest  level. 
At  the  refueling  contact  position,  two  main  vortex  effect  remains:  strong  downwash  and  pitching  vortex 
gradient. 


VI.  Conclusion 

In  this  paper,  the  equations  of  motion  of  a  receiver  aircraft  are  written  in  terms  of  its  position  and 
orientation  relative  to  the  tanker  and  its  inertial  velocity  as  the  sum  of  its  velocity  relative  to  the  surrounding 
air  and  the  velocity  of  the  air  relative  to  an  inertial  frame.  This  enables  the  formulation  of  the  relative  motion 
of  the  receiver  and  facilitates  the  incorporation  of  the  vortex  induced  wind  effect  into  the  translational  and 
rotational  equations  of  motion.  Engine  dynamics  has  been  modeled  as  a  first-order  system.  Saturation  and 
rate  limits  are  included  in  the  model  of  actuator  dynamics.  An  averaging  technique  has  been  outlined  to 
incorporate  the  effects  of  the  vortex-induced  wind  on  the  receiver’s  dynamics  using  the  standard  equations 
of  motion.  By  this  averaging  technique,  the  non-uniform  wind  distribution  induced  by  the  vortex-field  is 
approximated  as  uniform  effective  wind  components  and  gradients  acting  at  the  receiver’s  CM. 

A  multi-input  multi-output  position-tracking  state  feedback  controller  is  designed  so  that  receiver  aircraft 
tracks  its  reference  trajectories  specified  in  the  tanker’s  body  frame.  To  ensure  zero  steady-state  error,  the 
linearized  equations  of  motion  around  the  steady-level  refueling  flight  condition  are  augmented  with  three 
integral  states  for  the  position  tracking  error  in  the  three  axes  of  the  tankers  body  frame.  The  optimal  LQR 
method  is  applied  to  the  augmented  system  to  obtain  the  constant  gain  matrix  for  state  feedback  and  integral 
states.  The  controller  is  implemented  in  a  MATLAB/Simulink  based  simulation  environment  with  the  full 
6-DOF  nonlinear  equations  of  motion  for  both  the  tanker  and  receiver  including  the  wind  terms  as  well  as 
a  vortex  model  employing  an  averaging  technique  to  compute  the  effective  translational  and  rotational  wind 
velocities.  The  controller  is  not  provided  with  the  information  of  the  trailing  vortex. 

The  simulation  and  controller  are  applied  to  a  pair  of  tanker-receiver  aircraft:  KC-135R  and  ICE  aircraft 
(a  tailless  fighter  aircraft  with  innovative  control  effectors  and  thrust  vectoring  capability).  The  performance 
of  the  controller  is  evaluated  in  the  first  phase  of  a  refueling  operation:  the  receiver  aircraft  starts  away  from 
the  tanker  and  approaches  the  refueling  contact  position  by  tracking  a  commanded  trajectory.  Because  of 
the  redundancy  of  the  control  variables  (6  control  variables:  three  control  effectors,  thrust  magnitude  and 
two  thrust  vectoring  angles),  three  control  allocation  scenarios  are  simulated:  (1)  combination  of  control 
effectors  and  thrust  vectoring,  (2)  no  thrust  vectoring,  and  (3)  pitch  flaps  and  clamshells  are  fixed  at  their 
nominal  values  and  thrust  vectoring  used  for  pitching  and  yawing  moment  generation.  The  simulation  results 
show  that  the  controller  is  capable  of  tracking  the  commanded  trajectory  without  using  the  information  of 
the  nonuniform  wind  field  in  all  the  three  control  allocation  schemes.  However,  at  the  contact  position  the 
steady  state  condition  of  the  receiver  may  involve  non-zero  yawing  angle  and  in  the  third  case  non-zero  bank 
angle  due  to  the  vortex  induced  wind  components  and  gradients.  This  will,  in  the  next  phase  of  this  research 
work,  be  addressed  by  adding  two  more  integral  controls  for  the  yaw  and  bank  errors  since  there  are  enough 
number  of  control  variables.  The  simulation  is  also  proven  to  be  a  very  realistic  -  at  least  qualitatively  - 
evaluation  environment  to  determine  the  effect  of  the  trailing  wake  vortex  on  the  receiver  aircraft  and  to 
assess  the  performance  of  any  aerial  refueling  controller.  Another  item  of  future  work  is  to  investigate  the 
performance  of  the  controller  when  the  tanker  aircraft  is  in  a  ’’racetrack”  maneuver  when  the  tanker  starts 
to  enter  a  turn  with  a  max  bank  angle  of  30  degrees,  and  holds  this  position  until  it  is  180  degrees  from 
its  original  heading,  then  it  resumes  straight  flight.  The  receiver  must  be  able  to  remain  in  the  contact 
position  with  the  boom  connected  throughout  this  maneuver.  Moreover,  in  the  vortex  induced  wind  field, 
the  stochastic  content  of  the  wind  will  be  determined  and  modeled  so  that  the  performance  of  the  controllers 
can  be  evaluated  in  the  vortex  induced  wind  field  with  turbulence. 
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A.  State,  Control  and  Disturbance  Matrices 


The  state  matrix  of  the  linearized  equations  of  motion  at  the  steady-level  flight  condition: 
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The  control  matrix  of  the  linearized  equations  of  motion  at  the  steady-level  flight  condition: 
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The  disturbance  matrix  of  the  linearized  equations  of  motion  at  the  steady-level  flight  condition: 
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